Electromagnetic radiation from vortex flow in type-II superconductors 
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We show that a moving vortex lattice, as it comes to a crystal edge, radiates into a free space the harmonics of 
the washboard frequency, ujo = 2'Kv/a, up to a superconducting gap, A/ft. Here v is the velocity of the vortex 
lattice and a is the intervortex spacing. We compute radiation power and show that this effect can be used for 
generation of terahertz radiation and for characterization of moving vortex lattices. 

PACS numbers: 74.25.Qt, 74.25.Nf 

Effects at the washboard frequency in the flux-flow regime 
of type II superconductors were first reported by Fiory HI. 
Imposing an rf current at the harmonics of tJo = 2'Kv/a on 
top of the dc transport current, Fiory observed steps in the I- 
V characteristics of superconducting aluminum films. Here v 
is the vortex velocity and a is the vortex lattice parameter in 
the direction of v. Similar observations were made in high- 
temperature superconductors 12), |3(]. Larkin and Ovchinnikov 
and Schmid and Hauger js]] have demonstrated that the ef- 
fect originates from the coherent action of defects as they are 
passed periodically at the washboard frequency by moving 
vortices. It is analogous to Shapiro steps in Josephson junc- 
tions in the presence of the ac current [6j. By analogy with 
the radiation from a Josephson junction [7, 8, 9J, the inverse 
effect of the electromagnetic radiation from a moving vortex 
lattice at the washboard frequency should be expected. Radi- 
ation produced by a single vortex or vortex bundle crossing 
the edge of a superconductor was discussed by Dolgov and 
Schopohl They used the transition radiation approach to 
calculate the energy and spectrum of the broad-band electro- 
magnetic pulse emitted by the bundle. In this Letter we show 
that the situation for a moving vortex lattice is completely dif- 
ferent. As vortices come to the surface periodically with the 
frequency v/a, they radiate into free space the harmonics of 
ll>o up to a frequency corresponding to the superconducting 
gap, A/?i. Unlike the rf coinponent in the transport current 
produced by defects ([llSSBEtlj the electromagnetic radia- 
tion from the surface does not disappear but becomes stronger 
when disorder weakens. It is generated by oscillating electric 
and magnetic fields of vortices near the surface and propagates 
into free space due to the continuity of tangential components 
of the fields at the surface. We study cases of ideal and dis- 
ordered vortex lattice in small and large crystals as compared 
to the radiation wavelength. Radiation power is derived and 
shown to be within experimental reach. For a large crystal 
the effect is proportional to the radiating surface, while for a 
small crystal it is quadratic on the surface, that is, superra- 
diant. The challenge for experiment is to satisfy the require- 
ments of sufficient speed and appreciable correlation length of 
translational order. If these conditions are satisfied, the effect 
can be used for generation of electromagnetic radiation well 
into the terahertz range. 



FIG. 1 : (Color online) A moving vortex lattice radiates electromag- 
netic waves into a free space. Arrows indicate directions of the ap- 
plied magnetic field Ho, transport current j, vortex lattice velocity v, 
and the Poynting vector P^. Vortices are shown by circles. 



We study electromagnetic radiation from the boundaries, 
X = 0, —Lx, of a type-II isotropic superconductor, located at 
—Lx < X < Q, into a free space at a; > and x < —L^, 
see Fig. 1. The applied magnetic field, Hq, along the z axis 
is assumed to satisfy Hd <C Hq ^ Hc2, with Hd and Hc2 
being the first and the second critical fields. The transport 
current along the y-axis results in the motion of the vortex 
lattice at a speed v along the x--axis. We consider large v 
dominated by the Lorentz force and vortex drag, so that the 
effect of the surface Meissner current and the back effect of 
the radiation on the lattice motion can be ignored. To obtain 
the radiation power we need to find electromagnetic fields, 
Bz and Ey, at the surfaces a: = and x = —L^- At v ^ c 
the fields inside the superconductor can be obtained from a 
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quasistatic approximation (r = x, y, z): 



V X V X B + — B 



(1) 
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^ j dz5[v^v^{t,z)]Q[-x^{t)]Q[L 



c[v X E] = -ae/at. 



(2) 



where is the unit vector along the z-axis, <d{x) is a unit 
step function, A is the London penetration length, and rn(t, z) 
is the position of the center of the n = (ri,,p)-th vortex line 
in the vortex lattice. Consider first a square lattice with the 
lattice spacing a ~ ^J^q/ B <C A. We shall assume that the 
transport current j is sufficient to drive vortices in the laminar 
(Bragg glass) regime 111 Ull2ll : 



Xn{t) = an + vt , yp(z) = ap + Sp(z) . 



(3) 



Here n,p are integers and 5p{z) accounts for disorder in the 
position of vortices along the y-axis. Strictly speaking, vor- 
tices do not form an ideal periodic lattice along the x-axis. 
Rather, in the absence of topological defects, moving vor- 
tices adjust to quenched disorder by periodically replacing 
each other at some favorable points, x„ (0, p, z), with a period 
T = a/v m. 



Xn{t,P, z) = x„(0,p, z) + [a:;„+i(0,p, z) - a;„(0,p, z)]/(t). 

Here f{t) is a periodic function with the period T, satisfying 
f[t) K, t/T inside the interval [0,T]. One can show, how- 
ever, that at a <C A Eq. (|3]l is sufficient to describe the motion 
of Bragg glass along the x-axis - the magnetic field at the 
boundary (that determines radiation) is only weakly affected 
by disorder in a;„(0,p, z). Phase shifts in otherwise perfectly 
periodic process are introduced only by time-dependent ther- 
mal fluctuations. This is similar to Josephson oscillations. 



see Ref. [13|. The similarity of the vortex-induced radiation 
to Josephson oscillations also follows from the fact that the 
washboard frequency satisfies the relation hw = 2eV, where 
V is the voltage drop for one period of the vortex lattice in the 
direction of the current. This is a consequence of the expres- 



sion for the dc electric field in that direction, Ey 



(dc) 



{v/c)B. 



Note that delta-functions in the right hand side of Eq. ([T]i cor- 
respond to point-like vortex cores. In reality the cores have 
radius ^, the superconducting correlation length, which will 
be accounted for in the following. Here we will ignore ther- 
mal fluctuations that lead to the broadening of radiation lines. 

For 3> A one can prove that the two surfaces, a; = and 
at X = ~Lx, can be studied independently. We, therefore, 
begin with considering a semi-infinite superconductor at x < 
0. The fields inside the superconductor {x < 0), must be 
matched by solutions of Maxwell equations in the free space, 
{x > 0). The total fields inside the superconductor are 



Ey{t,r) = Eyy{t,r) + Eyo{t,r) , 



where B^v and Eyy are magnetic and electric fields produced 
by vortices, while B^a and Eyo are solutions of Eqs. ([T]i and 
(|2]i with a zero right-hand side in Eq. ([T]i. For Bzv one obtains 



= Y.fdz 



$0 



-iujan /v — ikyyp{z) 



771, p 



A^k^ iuj — ikxV — e 



«(1 + A2k2) 



iuj — ikxV — e 



(5) 



where k = {kx, ky,kz). The summation over n was carried 
out with the help of the relation 



M 



exp[-'{iu!a/v)n] = 27r 8{uja/v — 27rm), (6) 

n m— 1 

where the upper limit M ~ a/C is due to the nonzero size, ^, 
of the vortex core. The frequency dependent London length 
is given by A"^ = A"^ ^ ^5 + 47rfci^crg/c, where aq is the 
quasiparticle conductivity and k^ = oj /c. Integrating Eq. (|5]l 
over kx and using inequality a <C A we obtain the amplitude 
of the oscillating magnetic field at the boundary a; = at 
uj = Wm = muJo: 



B„,{u},x = 0,ky,kz) = 



p 



dz — — . (7) 

ziwAoj™ a 



(4) 



For the electric field of the vortex lattice, with the help 
of Eq. (|2]i and V • E = one obtains Ey,;{uj,'k) = 
{vk"^/ {c]i^)B2y{u!,'k). Solutions of the homogeneous equa- 
tions are Eyo{uj, x, y) — —ik^X^B zo{oJ , x, y) and 

BA^^x^y)^ j ^Aiu;,ky)e<^'-^'+''l^'''^+^''yy, (8) 

where A{lu, ky) is determined by the continuity of the tangen- 
tial components of the fields at the boundaries of the super- 
conductor 

Maxwell equations determine the relation between B^ and 
Ey everywhere in free space a; > 0. To find this relation 
at the boundary, x +0, we follow the derivation outlined 
in Ref. We assume that at a; > there is only outgo- 
ing electromagnetic wave. Then the electric field Ey{uj, r) = 
J dt exp{iujt)Ey{t, r) in free space is completely determined 
by its value at the boundary through [kj^ = {ky, kz)] 

Ey{u;,r)^ J ^Ey{i.,x = 0,k^)e^'^-'^ 

X {exp[*(fc2 _ kiy/hign{u;)x]eikl - ki) 

+ exph(ki - kiy/^x)]e{kl - kl)] . (9) 

The first term in the right hand side is due to radiated waves, 
while the second term comes from the waves decaying expo- 
nentially away from the boundary. Using this relation, Eq. (|2]i 
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and V • E = 0, we express the magnetic field in free space 
via Ey{u!, X ~ Q,y, z) and finally obtain the relation between 
Bz{uj, X = Q,y, z) and Ey{uj, x = 0,y, z) at the boundary: 

B,(u;,0,ki) = C(tc',ki)S;,(cc;,0,kx), (10) 
|A;„|e(fc2 _ki) *fc^e(ki-fc2) 



For the left boundary, x = —L^, one should reverse the sign 
of C in Eq. 

The above relations allow one to determine the fields B^o 
and E'yo and, finally, the total fields and Ey at the bound- 
ary as well as the Poynting vector of the radiation. Radiation 
power at frequency lo outside the superconductor at a; > (to 
the right) is given by 



nad(^) 



47r 
dkj. 
Try (2^ 



j dydzRc[Eye-"^*] Rc[B^e- 



(11) 



where integration is over I kx I < {k^^l. Using Eqs. f4\ and (fTOli 
we obtain for the amplitude of the oscillating magnetic field 



Eyiuj,0,k^) 



Eyy{uj, O.ki^) +iXujk^Bzv{uj,0,'k±) 
1 + iXujk^Ci'^, 



(12) 

with Xujk^ <C 1 and Eyy{Lu,0, kj_) < X^k^B^vi^J-, 0,kj^) at 
uj = muiQ < A/h. The power of the radiation into free space 
to the right or to the left of the superconductor becomes 



lad 



ck^ f {dk^/An'^ 



■.{\X^k^B,,iLO,0,k^)f 



8n J VfcS^ 
±2lm[X^k^B,,{Lj,0,k^)Ey,{uj,0,ki^)]}. 



(13) 



We see that the radiation power is slightly stronger to the right 
(in the direction of the vortex motion), though the difference 
is small. In the following we neglect this difference. 

Substituting the amplitude of the magnetic field at the 
boundary, B^vi^^ 0, k^) of Eq. Q, into the expression (fTJI i 
for the radiation power, and averaging over disorder in the 
vortex lattice, we obtain the power. 



dk 



± 



S{k 



±) 



(14) 



327rca2 J (27r)2 ^fcj^ 

S{k±) = ^ f —exp[ik^z - ikyyp{z)] , (15) 
p 

where 5'(k^) is the structural factor of the moving vortex lat- 
tice. It is peaked at fcz = and ky = 2TTn/a, and has the 
following properties: 5(0,0) = LyLz/a? ^ J dk±S{k±) = 
(27r/a)2. If translational correlations in the vortex lattice de- 
cay exponentially, then at Ly \ <C fcy.z ^ {o.ly,z)~'^ one 
should use S{ky,kz) = (y£z, where £y and £z are dimension- 
less transversal and longitudinal correlation lengths expressed 
in units of the intervortex spacing a (Iz ^ ^y) nlSll . 



Consider first the case of weak disorder, when correlation 
lengths exceed dimensions of the crystal, a£y ^ Ly and 
alz 3> Lz, and when the wavelength of the radiation is small 
compared to the crystal, c/cj^ <C Ly,L^. With account of 
the above properties of the structural factor, integration over 
kj^ in Eq. (fT4] l gives for the radiation power at = u)„i 



'Prad('J-'m) — LyL 



v'B' 
327rc 



(16) 



where we neglected insignificant dynamical terms in A^^. An 
amusing feature of the above result is that each electromag- 
netic mode up to cOm ~ A/ft satisfying conditions c/w„i <C 
Ly, Lz makes equal contribution, v'^B^ /{32Trc), to the radia- 
tion power per unit square of the radiating surface. The total 

power, summed up over all modes, is of order 'P^Id*^ I ^^Ly = 
[l/(16V2^)](z;Vc)-B3/2/7y2 



(A//327r)(u2B2/c) - [\/{\QV^)]{v'' /c)B^/''H'^i\ For 
weak disorder but small crystal (or lower frequency LUm), 
Ly, Lz <C c/ujrn, we obtain 



(17) 



Quadratic dependence of the radiation power on Ly and Lz 
is due to superradiance, that is, coherent radiation by vortices 
positioned within the radiation wavelength. 

Stronger disorder results in smaller correlations lengths, 
aly ^ Ly, {c/ojm) and a£z ^ Lz, {c/ojm)- In this case the 
power of the radiation is reduced: 



^rad(t^m) 



a'^iyi^Lu'^ 



(18) 



The coherent radiation takes place within the surface area of 
size {lya){lza), while the total radiation power is a sum of in- 
dependent contributions from such areas. Note, that ty, lz are 
correlation lengths for the movine vortex lattice. Due to the 
effect of dynamical reordering 111, 12, 1/7], they must be 
greater than the ones for the static lattice. Thus measurements 
of the electromagnetic radiation from moving vortex lattices 
provide information on translational correlation lengths and 
their dependence on velocity. 

Static square lattices have been observed in high- 
temperature superconductors : Lai 83Sro.i7Cu04+5 in fields 
above 0.4 T ill and YBaaCugOy in fields above II T HI . 
Square vortex lattices have been also observed in borocarbides 
12011 . There is little information, however, on whether the rect- 
angular symmetry is preserved in a vortex lattice moving at a 
high speed. For a triangular vortex lattice Xnp in Eq. ^ must 
be replaced by .T„p = an + vt + a[l + {—l)P]/A. ThenEq. Q 
acquires an additional factor [1 + exp(i7rTO)] so that only even 
harmonics, 2'mujo, of the washboard frequency are present in 
the radiation; their intensity being 1/4 of that for a square lat- 
tice. For a vortex liquid the integrals must be dominated by 
correlations at small distances, that is by ky of order 1 /a, and 
the radiation becomes suppressed by a factor d^uJ^^jc? , i.e., it 
is practically absent. 
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Energy burst from a bundle of A = LxLy/a? vortices ization of moving lattices. 



crossing the edge of a superconductor has been discussed by 
Dolgov and Schopol [10]. They correctly noticed that the radi- 
ation spectrum from the vortex lattice consists of discrete lines 
and can be superradiant. According to Ref. [lOt the transient 
radiation of a single bundle, as in the case of a single vortex, 
peaks at (a/V3A)ijJo, which is approximately the inverse time 
needed for the vortex field of diameter A to cross the crystal 
edge. We found completely different spectral power of the 
electromagnetic radiation from a moving vortex lattice. 

To obtain significant radiation power one needs high vor- 
tex lattice velocity v. This velocity is limited by dissipation 
(heating). Pulse technique helps to diminish heating and reach 
higher velocities llisl Il9ll up to the critical velocity v* aris- 
ing from Larkin-Ovchinnikov instability 12111 . The latter oc- 



curs due to the decrease of the vortex drag coefficient with v 
as ri{v) — ?7(0)(1 + / v*"^)^^ . Such a behavior of is 
caused by the escape of quasiparticles from normal cores due 
to the effect of the electric field. It leads to a negative slope 
in the I-V characteristics at velocities above v* . The latter 
depends on the inelastic scattering of quas iparticles and de- 
creases on cooling. Experimental study |22'] of this instability 
in Ndi.ssCeo.isCuOx gives v* ~ ^ m/s at 8 K in the field 
120 mT. At such velocity we estimate the radiation power as 
2 /iW/cm^ at ^ A/?i in the field of 1 T for a square vortex 
lattice in the case of weak disorder, Eq. (fTSI l. Much higher 
value of the critical velocity, v* ~ 1.2 • lO'^ m/s, was obtained 
in YBa2Cu307_i- films at 72 K in the field of 1.8 T Q. At 
such speed and field the radiation power from a weakly disor- 
dered square lattice would be of order 0.1 W/cm^. 

To obtain the efficiency of the radiation source, we write 
the dissipation power as Vdis = {BLxLy/^Q){7]v'^ Lz/2), 
where B = $o/a^, (BLxLy/^o) is the number of vor- 
tices, and rjv^ /2 is dissipation of energy per unit length of 
the vortex line. For a large superconductor the efficiency, 
^ ^rad*V^dis, in the casc of weak disorder, is given by 
r — ^qB / {32TrcriLx)- To estimate rone can use the Bardeen- 
Stephen formula for the drag coefficient, ?/ = Hc2^o/{PnC^), 
with pn being the normal state resistivity. At ^ lOA and 
B ~ 0.1Hc2, r can reach lO^**. This means that at a moderate 
vortex lattice velocity of 20 m/s and a cooling rate of 1 W/cm^ 
the radiation power of 0. 1 mW/cm^ can be achieved. Note that 
small r justifies the approach in which the back effect of the 
radiation on the motion of vortices is ignored. 

In Conclusion, coherent electromagnetic radiation should 
accompany the flux-flow state of a moving vortex lattice. The 
spectrum of the radiation has discrete character and extends 
up to the frequency corresponding to the superconducting gap. 
Thus, in principle, this effect can be used to generate radiation 
in the teraherts frequency range. Radiation power depends 
strongly on the vortex lattice symmetry, velocity and degree 
of translational order, providing a possible tool for character- 
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